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Abstract 

The renormalization group (RG) method is extended for gfobal asymptotic 
analysis of discrete systems. We show that the RG equation in the discretized 
form leads to difference equations corresponding to the Stuart-Landau or 
Ginzburg-Landau equations. We propose a discretization scheme which leads 
to a faithful discretization of the reduced dynamics of the original differential 
equations. 
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It is a fundamental problem in physics, especially, in statistical physics since Boltzmann, 
to reduce many degrees of freedom of a dynamical system to fewer degrees of freedom pre- 
serving the essential nature of the system. The reduced dynamics is often described by 
a few collective variables which represent slow and long-wave length motion of the system. 
There are several methods for the dynamical reduction such as the multiple-scale meth- 
ods (including the reductive perturbation method the average methods (including 
Whitham method [§]), the method of normal forms and so on. We are usually interested 
in asymptotic behavior of the system after a long time. Thus the problem is to obtain a 
dynamical reduction in asymptotic and global domains. 

Recently it has been recognized and emphasized by Chen, Goldenfeld and Oono that 
the renormalization group equation first developed in quantum field theory is a powerful 
tool for global and asymptotic analysis: They applied the RG equation a la Gell-Mann-Low 
to ordinary and partial differential equations and showed that the RG equation nicely gives 
a reduction of the dynamics describing slow and long-wave length motions of the system; the 
reduced dynamics is described by so called the amplitude equations. Afterwards the reason 
of the powerfulness of the RG equations was accounted for in the context of classical theory 
of envelopes More recently, it has been shown that the RG method can also naturally 
leads to the phase equations as well as the amplitude equations 

The purpose of this paper is twofold: One is to show that the RG method can be extended 
to discrete systems and leads to dynamical reduction of discrete dynamical systems (maps). 
A discrete system is described by a difference equation. We notice that the extension is highly 
non-trivial because the applicability of the RG method for differential equations is essentially 
relied on the local nature of differentiation P). Another purpose is to propose a discretization 
scheme for differential equations: It is needless to say that different discretization schemes 
of non-linear differential equations leads to different dynamical systems. We shall show that 
such a discretization-scheme dependence also exists for the reducibility of the dynamics and 
give a good discretization scheme based on the notion of the reducibility of the dynamics. 

Now let us take the following discrete system as a typical example of nonlinear discrete 
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systems, 

Xn+2 -2cOs9Xn+l + Xn = e/ (x^, ) , (l) 

where ^ is a constant and the function f{xn, Xn+i) contains nonhnear terms. We remark that 
the equation can be converted to a vector equation, X„ = F(X„) with X„ = ?/„ = Xn+i) 
and F(X„) = *(?/„, 2 cos — Xn + ef{xn,yn))- We have taken an example of second-order 
equation here, but the following discussion is applicable for higher-order equations. We 
notice that the unperturbed equation (e = 0) has a neutrally stable solution exp{±i9). ||TT[| 
Assuming that e is small, we shall apply the perturbation theory. To make the discussion 
definite, let us take a concrete form for /(x„, as; Xn+i) = aiXn + 0'2x\ + ^z^^- If 
we put 02 = 0, the resultant equation may be considered as a discrete version of the damped 
Duffin equation. 

Expanding x„ as x„ = + ex^^^ + e^xi^^ + ■ ■ let us try to find a solution which is 
valid around n = uq where riQ is arbitrary: x^*^ {i = 0,1,...) satisfies Lx^^^^ = 0, I/X^^^ = 
O'lXn^ + 020;^°''^ + a^x^^^^ , and so on. Here L = — 2 cos 9 ■ E + 1 with E is the forwarding 
operator, i.e., Exn = Xn+i- The lowest order solution may be written as 

4°) = A(no)e'"^ + c.c. (2) 

Here we have made it explicit that A may be a function of no; its functional dependence is 
yet unknown and will be determined by the RG equation, which determination constitutes 
the central part of the RG method. The first order correction x^^^ is given by 

x« = -ir^in - no)ia,A + 3a3\A\''A)e'^'"'^' + c.c, (3) 

where we have omitted non-singular terms which are proportional to exp{ikn9) only with 
k = 0,2, 3. Notice that there appeared a secular term proportional to n — riQ. We remark that 
the solution of the first-order equation is not unique; one could add any term proportional 
to xl^\ We have chosen the above form so that only new independent terms appear and the 
secular term vanishes ai n = uq, which assures the lowest order approximation is as good as 
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possible. Up to 0(e^),we have an approximate solution Xn = x^^ + exj^^^ = Xn{{A{no);no), 
which is only locally valid around n = Vrio; the validity as an approximate solution is lost 
due to the secular term as \n — uqI becomes large. It means that the naive perturbation 
expansion breaks, which is a well known fact. 

One may say, however, that we have a family of discrete curves Xn{A{nQ),nQ) with hq 
being a parameter characterizing the curves [Q. Each curve will become a good approx- 
imation for the exact solution around n = no if A{nQ) is suitably chosen. Therefore the 
"envelope" of the family of the curves may give a good approximation of the exact solution 
in a global domain. [|12| The "envelope" of x„((y4(no); uq) is constructed as follows: We first 
impose that 

AnoXn{A{no);no) = 0, (4) 

where A„q denotes the difference w.r.t. no. This is the basic equation of our method 
and we call it the RG/E equation. This is an equation to give no and A{no) as functions 
of no(n). The "envelope" x^ is given by inserting this solution to x„(A(no); no); x^ = 
Xn{A{no{n))]nQ{n)). However, since we are constructing the "envelope" which contacts 
with the local solutions at n = no so that the "envelope" give a good approximation, no 
should be n,i.e., uq = n: Notice that this choice nicely eliminates the secular term from x^. 
Conversely speaking, A{n) can be determined so that Eq.(4) gives no = n; this possibility 
is related with the "renormalizability" of the equation 0. Thus we have 

A^(n) = + 3a3|A(n)|2A(n)), (5) 



where A denotes the difference operator w.r.t. n [jT3[. This is the amplitude equation in a 
discrete form corresponding to Stuart-Landau equation for continuum systems, is 
thus given by 

xl = (A(n)e"' + c.c) + 2a2e| A(n) p + h.h.. (6) 
here h.h. denotes the higher harmonics. 



A significant point is that this function gives an approximate but uniformly vahd solution 
in a global domain. Let us show this in a general setting. Let X„ = X2„, • • • , Xdn) 

and F(X„, n) = *(Fi(X„, n), F2(X„, n), • • • , Fd(X„, n)); we assume that F(X„, n) is analytic 
with respect to X„. We consider the difference equation 

AX„ = F(X„,n). (7) 

We suppose that X„(W(no),no) is an approximate solution to the equation up to 
O(e^), where the d-dimensional vector W(no) denotes the initial values assigned at 
the initial time n = no- Here notice that no is arbitrary. The envelope function 
is given by X^ = X„(W(n),n), where W(n) is the solution to the RG/E equa- 
tion, AjiQX„(W(no), TT-o) = O(e^). Then one can show that X^ satisfies the orig- 
inal equation uniformly up to O(e^) as follows: AX^ = AX„(W(no), no)|no=n+i + 
A„„X„(W(no),no)|no=n = F(X„(W(no, no), n)|„o=n+i + C>(ef) = F(X„(W(n),n),n) - 
dF/dXri ■ A„oX„(W(no),no)Uo=n + 0(e"^^P) = F(Xf ,n) + 0{eP). Here we have utilized 
the RG/E equation in the second and last equality, and the analyticity of F in the third 
equality. This concludes the proof. 

We notice that the amplitude equation Eq.(5) is a first-order equation and a dynam- 
ical reduction is achieved in comparison with the original equation. In fact, in the polar 
representation A{n) — i?„exp(i(^„), i?„ and satisfy the following equations, 

Rn+l — Rn — ^/'^ ■ {0'\Rn + SOsi?^), (8) 

and </7n+i — fn — e/2 • cot 9{ai + SasRD, respectively. Notice that Rn is determined by the 
first order equation independently of (pn, which in turn is given in terms of Thus one 
may say that the second-order dynamical system Eq.(l) is reduced to a first-order one. As is 
seen from the derivation, this reduced equation has a universal nature as has Stuart-Landau 
equation. 

The first order equation for has simple qualitative properties depending on the signs 
and values of ai and a^. For example, when ai < but > 0, the equation is converted 
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to fn+1 = fn + afn{l - fn) with /„ = ^^803/ 1 Oi | i?„ and a = e\ai\/2. For < a < 1, the 
equation has a fixed point /* = 1, while for 1 < a < 1.246, the map shows a two-period 
behavior, and after that the map rapidly shows multiple-period behavior then eventually 
becomes chaotic. 

As another example of ordinary difference equation, let us take the one which is derived 
as a discretization of the Rayleigh equation: x + x = ex{l — 1/3 • x^). We remark that the 
equation admit a limit cycle with the radius of 2 |jl5|. We take the following discretization; 
X —>■ {xn+i — 2x„ + x„_i)/At^, X — > {xn — Xn-i)/At. That is, the central difference for the 
second derivative and the backward difference for the first derivative. Thus we have 

Xn+i - 2 cos6'x„ + Xn-i = e(x„ - x„_i)(l - ^ ■ ^^2""^^ (9) 

where cos^^ = 1— At^/2 and e = eAt. We remark that this difference equation has a neutrally 
stable solution as the unperturbed one. This is not the case for other discretization schemes 
such as X — s> {xn+2 — 2x„+i + Xn)/ At^ and x {xn+i — Xn)/ At. We put 1/ At = OJ. 
In the first-order approximation with respect to e, we have 

Xn ^ no) = \A{no) + e^^^^^(l - u^\B{no)\'')B{no) ■ {n - n,)Y^' 

2i sm 6 

_l 2 ^(no)^exp(3mg) 

a'^ (exp(3^0) - exp{te)){exp{3te) - exp{-td)) '''''' ^ ^ 



where -8(^0) = {exp{i6) — l)A{no). Now the RG/E equation A„„a;„(y4(no); 
gives the amplitude equation AA{n) = eexp{—i9)/2ism9- {l — Lj'^\B{n)\'^)B{n), accordingly, 



= 0, 

no=n 

21 Ut^\l2\ 



A{n + 1) = A{n) + zA{n){l - \A{n)\^), (11) 

with z = eexp{—i6)/ cos6/2. If we take the polar representation A{n) = Rnexp{i(fn), 
we have Rn+i = Rn + €a'Rn{l — R^), and (pn+i = V'n + esin^^/2 ■ (1 — R^) with a' = 
At cos 6 / {2 cos 9 /2) . We note that a' as a function of At is a parabola-like shape and takes 
the maximum about .6 at At ~ .9; it vanishes at At = and At ~ 1.4. Thus we see that 
Rn goes up monotonically to a fixed point 1. 

With this Rn and (fn, the envelope x^ is given as 



= XniA{n);n) = 2RnCos{n9 + (fn) + 



e tan^ 



Rlsm{3{ne + ipn)~^e}, (12) 



12cos3e/2 



which shows that the radius of the hmit cycle is 2, irrespective of the choice of At in 
accordance with the original Rayleigh equation. We remark that it is not the case for other 
discretization schemes as given below Eq.(9). 

Fig. la shows given by Eq.(12) and the envelope 2i?„ together with the exact solution 
of Eq.(9) with e = .4 and At = 0.25. One can see that the agreement is excellent in the 
global domain; notice that the results is obtained in the first order approximation. One 
can also see that the amplitude 2i?„ successfully describe the slow motion of the system. 
The characteristic features of the system as a dynamical system may be more clearly seen 
in Fig. lb where the behaviors in the "phase space" are shown. The agreement is excellent 
again. [|1^] 

Finally, we consider partial difference equations. As an example, we take the difference 
equation which is given by a discretization of the 1-dimensional Swift-Hohenberg equation 



]T7| ; dt(f){x, t) = e{(j){x, t)—(j){x, t)^) — (l+9^)^0(x, t). With the discretization dt(f) {4>{n, m+ 
1) — (f){n, m))/At = Am4>{n, m)/At, dlcj) — > (0(?7, + 1, m) — 20(r;,, m) + 0(n — 1, m))/Ax^ = 
A^0(n, m)/Ax^, we have the following partial difference equation 



equation admit a dynamical reduction giving an amplitude equation which is the analogue 
of the time-dependent Ginzburg-Landau equation in the continuum theory. The reason why 
the dynamical reduction is possible in the RG method is that the difference equation is 
constructed so that the unperturbed equation has a neutrally stable solution. 

Making the Taylor expansion = 0o+e0i+e^02 + " ■ ■, we have equations in the successive 
order £0o = 0, £0i = 0o — 00, and so on. We consider an asymptotic solution at m — > oo 
and take the following neutrally stable solution as the 0-th order one; 



C(f){n, m) = e{(f){n, m) — (f){n, m)^) 



(13) 



where C = A^ + r(Ax^ + A^)^, with r 



At/Ax^. We shall show that the difference 



(poin^m) = A (no, 7710)6*"^" +C.C, 



(14) 
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with 6x — 2Sin~^Ax/2. Here we have made it exphcit that the amphtude A may depend in 
a yet unknown way on the initial time mo and point uq. They will be determined by the 
RG/E equation. 

Then the first-order equation now reads 

£01 = {{A - 3|A|2A)e''^^- - A^e^*"^-} + c.c. (15) 

A straightforward but somewhat tedious manipulation gives 



)i(n, m) 



IXi{m - mo) - TTT^-TTrin'^^^ - no' + ie *^"(n - no)} 
8sm i)x 

x(A-3|A|2A)e''*^- - 



+C.C., (16) 



where rS"^^ = n{n — 1) and fii + fi2 = 1- Thus we have an approximate solution (f)(n,m) = 
4>o{n, m) + e4>i{n, m) which is valid only for (n, m) around n ~ no and m ~ mo- 

Now the RG/E equation /^mo(t>\mo=m = and A„o0|no=n = gives A^A(n, m) = e//i(l - 
3|A|2)A, and A,2y4(n,m) = (1 -3 jAp) A, respectively. Her c wc have utilized the fact 

^mA = 0(e) and A„74 = 0(e), and neglected terms of O(e^). Thus noting that ^1 + ^2 — ^1 
we reach the amplitude equation for the difference equation 

A^^(n, m) = 4r sin^ 6I^A^A(n, m) + e(l - 3|A(n, m)\^)A{n, m). (17) 

This is precisely the discretized form of the time-dependent Ginzburg-Landau equation: 
dtA{x, t) = Adl(t){x, t) + e(l - 3|A(x, t)\'^)A{x, t). 

In summary, we have shown that the renormalization group method can be nicely ex- 
tended to discrete systems, and that the method is useful as a tool for global asymptotic 
analysis and gives dynamical reduction of discrete systems. We have emphasized that the 
method is applicable for systems which have neutrally stable solution. It is to be remarked 
that this is also the case for equations which have unperturbed solutions on invariant sta- 
ble and unstable manifolds. We have also given a notion of the discretization scheme that 
faithfully preserves the nature of the reduced dynamics irrespective of the magnitude of At. 
Finally, we notice that it is not trivial what discretization scheme preserves the integrability 



of differential equations which admit sohton solutions [0. It will be interesting to examine 
if our discretization scheme based on the RG method can have a relevance to soliton theories. 
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Figure Captions 



Fig.la The dots shows (Eq.(12)) v.s. nAt while the thin hne the envelope 2i?„ for 
€ — A and At — 0.25. The bold hne shows the exact solution of Eq.(9). 
Fig. lb The behaviors in the "phase space"; the vertical axis denotes the "velocity" Vn = 
{xn+i — Xn)/At, while the horizontal QjXIS Xfi . The dots shows our approximate solution 
whereas the solid line the exact one. 
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Fig. la 
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